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Abstract 

We study lepton flavor violation (LFV) in tau decays induced by heavy 
Majorana neutrinos within two models: (I) the Standard Model with addi- 
tional right-handed heavy Majorana neutrinos, i.e., a typical seesaw-type 
model; (II) the Standard Model with left-handed and right-handed neu- 
tral singlets, which are inspired by certain scenarios of S'O(IO) models and 
heterotic superstring models with E§ symmetry. We calculate various LFV 
branching ratios and a T-odd asymmetry. The seesaw Model I predicts very 
small branching ratios for LFV processes in most of the parameter space, al- 
though in a very restricted parameter region it can reach maximal branching 
ratios Br{r — > ^7)~10 -9 and Br(r — ► 3^)~10 -10 . In contrast, Model II may 
show branching ratios Brir — > e7)~10~ 8 and Br(r — > 3e) ~ 10~ 9 , over a siz- 
able region of the parameter space, large enough to be tested by experiments 
in the near future. 
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I. INTRODUCTION 



One of the many puzzles remaining in the current phenomenology of particle physics 
is to understand the smallness of the masses (~ 1 eV) of standard neutrinos u e , and 
u T , compared to those of charged leptons. If neutrinos are of a Dirac nature, nonzero 
masses could be obtained in the Standard Model (SM) by introduction of (sterile) right- 
handed neutrinos. On the other hand, if neutrinos are of a Majorana nature, more appealing 
solutions to the small neutrino mass problem exist. In order to avoid the explicit breaking 
of the SM gauge symmetry and still obtain nonzero Majorana mass terms (via spontaneous 
symmetry breaking), an additional Higgs triplet is needed in the SM. The latter would result 
in physical Goldstone bosons, but these have been excluded by experiments at LEP. On the 
other hand, various models in the context of extended gauge structures result in Majorana 
mass terms and could give possible solutions to the neutrino mass problem. An appealing 
solution is the seesaw mechanism |l| within the framework of 5*0(10) or left-right symmetric 
models. In the conventional seesaw models, the effective light neutrino masses are within 
the scales of eV to MeV via a relation involving the hierarchy between very large Majorana 
masses and Dirac masses comparable to those of charged leptons. Another possible solution 
was investigated in the framework of heterotic superstring models with E& symmetry or 
certain scenarios of SO(10) models ||, where the low-energy effective theories include new 
left-handed and right-handed neutral isosinglets and assume conservation of total lepton 
number in the Yukawa sector. 

One possibility to test the neutrino sector lies in the study and measurement of lepton- 
flavor-violating (LFV) processes, e.g., /i — > or 3e; r — > 717 or 3/i; r — > or 3e. Such 
processes are practically suppressed to zero in the SM, due to the unitarity of the lep- 
tonic analog of the CKM mixing matrix and the near masslessness of the three neutrinos. 
Motivated by the aforementioned models with an extended neutrino sector, the authors 
in Refs. derived analytic expressions for LFV decay rates of charged leptons in such 

models with heavy Majorana neutrinos. The authors of Ref. |F] gave a model-independent 
framework for analyzing /i — > e'j and fi — > 3e processes and investigated specific features 
of several Supersymmetric GUTs. They focused on Parity- and T-violating asymmetries 
involving muon polarization in the initial state. 

Some generic properties of LFV processes and the corresponding constraints on the 
neutrino mass matrix have been studied in Ref. H. Phenomenological studies of various 
LFV and lepton-number-violating processes have appeared in the literature, including direct 
production of heavy Majorana neutrinos at various colliders 0, heavy Majorana mediated 
processes flOH , and LFV processes (including [i — > ej and r — > /ry) in supersymmetric 
frameworks [ IT ]. 



In this paper we will consider LFV decays of tau leptons in the framework of the two 
aforementioned models with extended neutrino sectors. We will concentrate on the calcu- 
lation of the corresponding LFV branching ratios and the corresponding expected numbers 
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of events. In addition, we will calculate a T-odd asymmetry induced by these processes.^] 
In Sec. II we review the two models in question. In Sec. Ill we present the formulas for the 
branching ratios of charged lepton decays, I — > l >r y and I — > 31', and the T-odd asymmetry for 
/ — > 31' within these two models. In Sec. IV we present approximate maximal values for LFV 
tau decay rates, exploring the possibility of obtaining sizable rates that can be measured in 
the foreseeable future, yet keeping consistency with present experimental constraints. We 
give a summary and state our conclusions in Sec. V. 



II. TWO NEUTRINO MIXING MODELS 



To set up our notation, we briefly review the two models in question: We call Model I the 
SM with the addition of right-handed neutrinos (singlets under the gauge group) and with 
the seesaw mechanism involved, and Model II the SM with the addition of both left-handed 
and right-handed neutral singlets. 

Model I: It is the SM with its Nl standard left-handed neutrinos vl% and an additional set 
of Nr right-handed neutrinos ujh, where the neutrino mass terms (after gauge symmetry 
breaking), which can be written as 

-C^=\{vL^M^ L \+h.c. , (1) 

contain a (Nl + Nr) x (Nl + iV^-dimensional matrix M. with a seesaw block form [|TJ. 
This matrix can always be diagonalized by means of a congruent transformation involving 
a unitary matrix U, namely: 

M = ( ° T mD ) , UMU T = M d . (2) 

The resulting N^+Nr mass eigenstates rij are Majorana neutrinos, related to the interaction 
eigenstates v a by the matrix U 

/ . . \ N L +N R / c \ N L +N R 

c = E Kn Li , [ Vl ) = E U ia n Ri . (3) 

\ R / a i=l \ R J a i=l 

The first Nl mass eigenstates are the light standard partners of the charged leptons, while 
the other Nr eigenstates are heavy. It is convenient, as done in Ref. ||, to introduce a 



1 Under the assumption of CPT symmetry, CP violation is equivalent to T violation. While 
standard CP violation appears in the quark sector, it could also arise, for example, in processes 



which involve only elementary (SM) bosons or (heavy) Majorana neutrinos 
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Nl x (Nl+Nr) -dimensional matrix B for charged current interactions, and a (Nl+Nr) x 
(Nl + Nr) -dimensional matrix C for neutral current interactions 

N L 

B h = U*, C l3 = Y,U ia U* a) (4) 

0=1 



where the charged leptons are taken in their mass basis. The ratio between the Dirac mass 
{mo) and the Majorana mass (tum) characterizes the strength of the heavy-to-light neutrino 
mixings (s^) 2 = J2h \Um\ 2 (~ I^dIVI^m] 2 ), as well as the size of the physical light neutrino 
masses: rn Ulight ~m 2 D /mM- In this model the very low experimental bounds on m vu ht (~ 1 
eV) impose severe constraints on the \mjj\ <C \ttim\ hierarchy required, and consequently 
also on the heavy-to-light neutrino mixings. 

Model II: It is similar to Model I, except that it contains an equal number Nr of left- 
handed (Su) and right-handed (vj&) neutral singlets 0H, and the form of the mass matrix 
M. is such that total lepton number is conserved (although lepton flavor mixing is still 
possible). After electroweak symmetry breaking, the neutrino mass terms are 

i I m D 

-C^r = -(u L ,V%,S L )M u R \+h.c, M=\m T D m T M \. (5) 
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rriM 



The mass matrix A4 is (Nl + Nr) x (Nl + Nr) -dimensional, where Nr = 2Nr (the Dirac 
block rri£) is A^^ x iV^-dimensional). When Nr = Nl, this model predicts, for each of the 
Nl generations, a massless Weyl neutrino and two degenerate neutral Majorana neutrinos. 
Consequently, the seesaw-type restriction m UUght ~ m 2 D /m M of Model I does not apply 
here p!5|JIE| . Here it is not the smallness of light neutrino masses but the present experimental 
bounds on the heavy-to-light mixing parameters (s u L 1 ) 2 ~ |^z)| 2 /l m Af | 2 (~ 10~ 2 , see below) 
which impose a certain level of hierarchy \uid\ < \ttim\ between the Dirac and Majorana 
mass sector. This hierarchy is in general much weaker than in seesaw models. Although 
Model II features (Nl) massless neutrinos in the light sector, nonzero masses for the light 
neutrinos can be generated by introducing small perturbations in the lower right block of 
A4, i.e., small Majorana mass terms for the neutral singlets Sui without much effect on the 
mixings of heavy-to-light fields. 



III. FLAVOR VIOLATING TAU DECAYS WITHIN THE TWO MODELS 



Recently LFV processes have been investigated extensively because SUSY GUTs predict 
that the branching ratios for fi — > ej and \i — > 3e and the fi — e conversion rate in a nucleus 
can reach just below present experimental bounds (7,17]. Here we address the predictions 



for LFV decays of the form I — > I'^y and I — > 31' within the models of Section II. 

The amplitudes for I —>■ l >r y and / — > 3/' in terms of the model parameters were obtained in 
Ref. These processes occur only at one loop level or higher in the (extended) electroweak 
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theory. The amplitude for / — > i'7 arises from a 7-penguin with the photon on the mass 
shell and is given by an expression of the form: 

M(l^ I'j) = iT^r^ uvia^imyPL + m,P fl )«, , (6) 



while the amplitude for I — > 31' receives contributions from 7— penguins, Z— penguins and 
box diagrams, M(l — > 3/') = M-y + Mz + Mbox- 



M^l - 


+ 3f) 


M z (l- 


+ 3Z') 


M B ox{l - 


+ 3Z') 



• ^"W w - 



2M 



w 



F l h lii - q -4)P L - G ll l lo ,X-{m v P L + m l PAu h (7) 



a. 



8M 2 , 



u v1v PlUx uvf (2 - As 2 w )P L - As 2 w P R v v F$ 



(8) 
(9) 



In the above expressions, e M [Eq. (H)] is the photon polarization, Pr/l = (1 ±7s)/2, and the 
factors G 1 ^ , -Fb 0;c are combinations of mixing matrix elements and some special functions 
that appear in the loop diagrams of the corresponding processes: 



Gt 



M' 



pll' 
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Box 



N L +N R 
H B li B l'i C 7 (Aj) , 

8=1 

N L +N R 

B uBi'i F-y(Xi) , 

i=i 

N L +N R 

B* u Bvj 5ij Fz{K) + Cij H z {Xi, Xj) + C*- G z {Xi, Xj) 

i,j=l 
N L +N R 

[^Bi'iBi'jB^Bi^ FBox(Xi, Xj) + Bi'iBi'iBijBi'/j Gsox(Xi, Xj 

ij=l 



(10) 

(11) 
(12) 

(13) 



The explicit expressions for the loop functions G 7 (x), F Box (x,y) are given in Ref. [|J, 
and their arguments are A, = mj/M^, i.e. the masses (squared) of the Majorana neutrinos 
inside the loop, in units of Mjy. Eqs. (|iT)|)-(|T3"|) involve a summation over all Majorana 
neutrinos, N R being the number of heavy ones (=N R , 2N r in Models I, II, respectively). 

From the amplitudes (|6])-(|9]), the decay rates are obtained. Using the notation of Ref. 0, 
they take the form 

T(l -> r 7 ) = T(l -> eP e Pi) ■ 38An 2 (l +ml/mf)\A R \ 2 , (14) 
T(l -> 3/') = r(/ -> eP e ^) • |2|(7 4 | 2 + |^ 6 | 2 + 8i?e (eA fl (2(/2+<£)) 

+ (l+m 2 ,/m 2 ) (321og (m 2 /3m 2 ,) - 104/3) \eA R \ 2 \ x jl + £>(m 2 ,/m 2 ) j , (15) 
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where g4,g& and Ar are the coefficients of the operators in the effective lagrangian relevant 
to these processes, and are given by: 

9, = ^{2s 2 w if + (f - 24) if + F*J , (16) 
g 6 = ^{2sl if + (-24) if} , (17) 

e ^ = l£ < > (18) 



where = g|/(47r) = V^G^M^/tt and a em = e 2 /(47r) are the weak and electromagnetic 
fine structure constants, and 4 = sin 2 ^. 

A T-odd asymmetry can be defined in the decays I — > 31' which is sensitive to the 
CP phases of the neutrino mixing matrices, but has the experimental drawback that it 
requires independent knowledge of the initial lepton polarization (in this case, the tau lepton 
polarization). Defining, in the CM frame, the decay plane as the plane that contains the 
three final momenta, At is the asymmetry between the cases when the polarization of the 
initial lepton points to one or to the other side of the decay plane. Geometrically, Ax is 
a </>-angle asymmetry, where <ft is the angle between the decay plane and the plane that 
contains the polarization vector of the initial lepton I and the momentum of the final lepton 
with charge opposite to / (see Ref. for details). The explicit expression for At is then: 




and in terms of parameters (|T6| ) -(|T8| ) it is 

At = V TuSm \j^ Im {eAR 9l) - ii Im {eAR 9l) } x I 1 + °^'i m ^\ ( 20 ) 



IV. NUMERICAL RESULTS AND DISCUSSIONS 



Several experiments provide constraints for the masses and mixings of light and heavy 
neutrinos: Tritium beta decay provides the present bound on the electron neutrino mass 
m Ve < 3 eV 1 18]. The solar neutrino deficit [19] can be interpreted either by matter en- 
hanced neutrino oscillations if Am^ ~ 1 x 1CT 5 eV 2 with small or large mixing, or by 
vacuum oscillations if Am 2 o/ ~ 10~ 10 eV 2 with maximal mixing ||20|1 . Atmospheric neutrino 



experiments show evidence for Am 2 im ~2.2 x 10 3 eV 2 with maximal mixing pl| , p2| . We 



will assume that Am 2 oi 



2 2 

\mt —rat 



and Am 2 tm 



I 2 2 

Imf, —m„ 



Am 2 tm as well. Since Am 2 



then 

to all three light neutrino masses: 



\ m v T ~ 
)2 



m; 



Since Am 2 oZ « Am 2 tm . 



m Ve ,m 



atm 

m,. 



<<3 2 eV , the 3 eV upper bound applies 
< 3 eV. Experimental evidence indicates 



that v T -v e mixing is (nearly) zero [ 23| . Further, Refs. |23| investigated possible patterns 
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of the Majorana neutrino mass matrix which are compatible with these results and the 
non-observation of neutrino less double beta decay [24]. In the models we are considering, a 



number of low-energy experiments set upper bounds on possible non-SM couplings, which 

[j4l,[26f as (s u L 1 ) 2 = J2h \Bih\ 2 (where h indicates heavy neutrinos). 



are characterized in Refs. 
Recent analyses |[27||, for models where the additional neutrinos are S{7(2)i-singlets, give 



< 0.005, 



{s^f < 0.002, 



< 0.010. 



(21) 



There is also a theoretical constraint, a perturbative unitarity condition (PUB) ||28|| , which 
states that perturbation theory to one loop is applicable only if the decay width Y nh of a 
heavy Majorana neutrino is small compared to its mass, say, T nh < M nh /2. In the limit 
of large masses M nh » Mw, Mz, Mh, the PUB constitutes an upper bound for heavy 
neutrino masses: ||,29,30| 



M: 



N L 
1=1 



Blh\ < , h = l, 

OLid 



N 



R- 



(22) 



In addition, there is a lower bound, M nh > 100 GeV, arising from the non-observation of 
heavy neutrinos in experiments to date. 

As the bound on (s^) 2 in (PH) is tighter than those on (s£ e ) 2 and (s^ T ) 2 , LFV muon 
decays are more suppressed than tau decays. We will therefore study LFV tau decays in the 
two models, trying to see if the present experimental upper bounds on r — > Ij and r — > 31 

Br{j- -> e~i) < 2.7 x 10" 6 , Br{j- -> pTi) < 1.1 x 10" 6 , (23) 
Br{r- -> e~e~e + ) < 2.9 x 10" 6 , Bt{t~ -> \i~ \T < 1.9 x 10" 6 , (24) 



can be reached theoretically once we account for all the aforementioned constraints. 

A numerical analysis of such reactions in Model I, for Nl = 3 and Nr = 2, has been 
performed in Ref. Q[. A comprehensive numerical analysis in Model II has been performed 
in Ref. J}1]]. In these two references, the analyses were performed by starting with the 
neutrino eigenmasses and specific combinations of the mixing matrix coefficients B^ on 
which restrictions were imposed. Our approach will be somewhat different, starting with 
explicit mass matrices, Eqs. @ and (f|), and from there deriving the masses and mixings. 
This approach is cumbersome if we include all three light generations, so we will take Nl = 2 
and Nr = 2. In this way, Brir — ► /i/ie, ee/i) will not be considered. However, since we are 
particularly interested in the largest possible branching ratios, and since either Brir — > 3e) 
or Br{r — > 3/x) is very suppressed (as argued below), then Br{r — ► fifie, ee/j) is also expected 
to be suppressed in comparison with the largest of the LFV branching ratios. 
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A. Model I 



In the considered case (Nl : 
generations (ui,£~) and (u t ,t~ 



phases in these types of models was studied in Ref. |L3 
two independent physical phases (Si, 



2 and Nr = 2) for LFV r decays, we have two light lepton 
where I is either e or \i. The structure of the CP-violating 

In the considered case there are 



1,2), so that the Dirac and Majorana submatrices 



in M. [see Eq. (EI)] can be taken to be of the form 

be iSl \ 



m D 



ce 



a 



d 



m M 



Mi 
M 2 



(25) 



where a,b,c,d are real. We take the convention M 2 > M\. The matrix M. can be diago- 
nalized via the congruent transformation of Eq. (0) - in numerical calculations we use the 
diagonalization approach as described in Ref. 
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We then find the values of mp (i.e. a, b, c, d; Si, 62) which give, for given heavy Majorana 
masses Mi and M 2 , the largest possible LFV branching ratios Br(r — > ji) and Br(r — > 3£). 

In Model I, the transformation matrix U of Eq. (|2|) can be presented as a product of 
a seesaw transformation block matrix U s and a light-sector mixing matrix V': U — VW S . 
The seesaw transformation U s produces an effective light neutrino mass matrix vn Vlight w 
m0Wjyinjj in the case mj) <mM, namely: 



m 



v light 



' a? 1 b 2 2i5i 



I — e 



+ 



M 2 C 



( ac_ i& 2 1 feti jgj 

VA/i C t m 2 c 
Mi c ^ Mi) 



) (i 



(i + o( 



m n m 



D" b M 



(26) 



and LFV mixings of order ~mom M . The light sector mixing matrix V*, which is the upper 
left part of Z7, is approximately unitary and of the form: 



Ft 



cos# sin^exp(— ze) 
sin^exp(ie) cos# 



(27) 



where ^ = and ir/4 correspond to zero and maximal mixing, respectively, and where e is 
a CP phase, in general a complicated function of Si and <$2- 

In the case r — > fi, we indeed have maximal mixing 9 — > 7r/4, according to atmospheric 
neutrino experiments. We will consider first this case. If we demand that this maximal 
mixing is obtained independently of the values Mi and M2 of the heavy Majorana sector 
[see Eq. (|25D1, then the following simple relations in the light Dirac sector are implied: 
a 2 = c 2 , b 2 = d 2 , and Si = S2 = S. The value of S can be restricted to lie in the range 
— Ti / 2 < S < 7r/2. The eigenmasses of the two light neutrinos are then 



m 



V\,Vl 




Mi M 2 



1/2 



± 



ac bd 
M~i + M~ 2 



(28) 
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while the heavy-to-light mixing parameters (|2~ID (s^) 2 = Ylt=3 \Bih\ 2 are 



(4-) 2 = W) 2 = ^ + -|(-l), (29) 



and the CP-violating parameter e of Eq. (|27|) is 



P a 2 /Mi - (b 2 M 2 ) , . 

tane = tan<5 x ) _ ' . 30 

a 2 /Mi + 6 2 /M 2 v ; 



Now, conditions a 2 = c 2 , b 2 = d 2 mean two possible cases: (1) a = ±c and b = ±d; or 
(2) a = ±c and 6 = =F<i. 

1. Case a = ±c and b = ±d: then > (a 2 /M 1 + b 2 /M 2 ), so s| < m v jM x < 3eV/M 1 < 
3 x 10 -11 . Since Br{r — > fi'-f) and i?r(r — > 3/i) are approximately proportional to 
( s l ) 2 (^l ) 2 (= ^l 4 ), ^ follows that Br(r -> /x 7 ) and 5r(r -> 3/x) are below 10" 24 . 
For muon decays, Br(fi — > e7) and I?r(/i — > 3e) are obtained by dividing the previous 
values by Br{r — > fiU^u T ) 0.174, thus obtaining 5r(// — > e7) and .Br(/i — > 3e) below 
10~ 23 , values which are well below their respective present experimental bounds (10~ n 
and 10" 12 ). 

2. Case a = ±c and b = =F<i: then m„ r > 2|a 2 /Mx — b 2 /M 2 \, with the equality being 
reached only when 5 = ir/2. In the latter case, m v = 0, and m Vr = 2\a 2 /Mi — 
b 2 /M 2 \ = (Am 2 ^) 1 / 2 « 0.047 eV. This case (<5 = 7r/2) thus avoids the suppression of 
s| = {a 2 /Ml + b 2 /M^) while keeping a 2 /Mi extremely close to b 2 /M 2 (a fine tuning 
situation). The value of s| can then be saturated to (s|) max — 0.002 [Eq. (plj)] with 
the following parameters in the Dirac matrix m^: 

a = c = M 1 {s L ) m J^l + MjM 2) b = -d = a(l±r l )^M 2 /M 1 , (31) 
MA 1 . „ -ii I f M l 



AmL 1 + - 2 « 1.17 x 10~ 13 x — f- I 1 + ^ ) , (32) 
V M 2 / 4M 1 (s 2 J max (s£) max V M 2 / 

and where, as mentioned, 5 = tt/2. The rates r — > /ry and r — > 3/i are again practically 
proportional to (sfj and approach their maximum (for fixed chosen values of Mi and 
M 2 ) in the case given by Eqs. (]3T|)-(|32|), as shown in the Appendix. The conditions 
(|3T])-([^), which are only reached by fine tuning, give the largest possible branching 
ratios in Model I, Br{r — > /i7)~10~ 9 and Br(r — > /i/i/i) ~ 10 -10 . In Fig. |l| we show 
the two branching ratios as functions of M 2 , for two different ratios M\/M 2 = 0.1 
and 0.5, accounting also for the PUB conditions (f22f). The CP-violating asymmetry 
parameter At (|20| ) is in this case, unfortunately, equal to zero, since 5 = tt/2 implies 
e = [Eq. (p0|)1 and thus no CP violation. If we move 5 away from 7r/2, the allowed 
branching ratios drop sharply, mainly due to the upper bounds m^, m UT < 3 eV, i.e., 
a situation similar to case 1 sets in. Accordingly, we do not consider other situations 
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of CP violation in Model I, as the branching ratios fall dramatically to unobservable 
values outside the fine-tuning condition. 

The results for Br(fi — > ej) and Br(fi — > 3e) are again obtained by dividing the above 
results by 0.174. However, we will then obtain values above the present experimental 
bounds 1.2 x 10~ n and 1.0 x 10~ 12 , respectively. We are thus led to conclude that the 
assumed fine-tuning condition is not really met in this case. 

If we now consider the case r — > e, i.e. the processes r — ► ej and r — * eee, the neutrino 
oscillation experiments indicate that the mixing is almost zero |23|]: 9 ~ in Eq. (p7|) . If 
we assume that the zero mixing condition is fulfilled independently of the heavy Majorana 
sector, we obtain the relations ac = and bd = 0. The cases where a = b = 0oic = d = 
give us (s£ e ) 2 = and (s v £) 2 = 0, respectively, and thus extremely suppressed branching ratios. 
The cases where a = d = 0oib = c = give: (s u L e ) 2 (s u L T ) 2 < m v lM x < 3eV/100GeV = 
3 x 10 -11 , i.e., as in case 1 discussed previously we obtain extremely suppressed branching 
ratios. 



B. Model II 



The neutrino mass matrix has the form of Eq. (|5]). In the considered two-generation 
case (Nl = 2, Nr = 4) for LFV r decays, the submatrices and can be taken in the 
form 

mo = cS d , m M = I Q Ma . (33) 



In the two-generation scheme of Model II there is only one CP-violating phase £ |fTo|. 

Since the (Nl) light neutrinos in the model are massless, the LFV r decay rates will 
neither be affected by the experimental light neutrino mass bounds, nor by the solar and 
atmospheric neutrino experiments and their requirements of the maximal {y IJi —p T , v e -v^) or 
minimal {y e —i> T ) mixing.^ However, the rates will be affected by the PUB restrictions 



as well as by the mixing parameter bounds (21) as the rates are proportional to {s v £) 2 {s v £) 2 
or (s v l) 2 (s v £) 2 . The r — > fi rates are suppressed in comparison to r — > e rates, because 
the upper bound for (s^ 1 ) 2 is smaller than that for (s£ e ) 2 [cf. Eq. (pip ]. Therefore, we will 
consider r — ► ej and r — >• 3e LFV rates. 



2 Nonetheless, it is possible to obtain nonzero light neutrino masses in Model II to accommodate 
neutrino oscillation experiments. Introduction of small mass terms for the neutral singlets Sn gives 
non-zero and non-degenerate masses of the light neutrinos, thus the possibility to accommodate 
Am 2 tm without significantly affecting the presented LFV rates. 
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Similarly as in Model I, we first find the Dirac parameters in Model II such that the LFV 
branching ratios, at fixed Majorana masses Mi and M 2 , are maximal. This occurs when 
inequality ( |A.3| ) in the Appendix becomes equality, and the values of the mixing parameters 
{sl) 2 are maximized, i.e., saturated according to Eq. (^l|). In the case of no CP violation 
(£ = 0), the requirement ( A.4 ) for the inequality ( |A.3| ) to become an equality gives the 



relation ad = be, while the saturation of the values of {s v £) 2 and {s v £) 2 gives two other 
conditions, for the four Dirac parameters a, b, c, d. This still allows us the freedom of fixing 
one of the four Dirac parameters without affecting the rates. We can, for example, require 
the symmetry of the (real) mo matrix: b = c. All of the above results in the following 
approximately "optimized" choice of mo parameters (when £ = 0): 

M 2 s lm (34) 



,2 _ „2 ' 
'lm °2m 



^/(M 2 /M 1 ) 2 + (s 2m /s lm ) 2 yfi 
b = c = a x (s 2m /s lm ) , d = a x (s 2m /s lm ) 2 , (35) 



where s\ m = (s^max = 0.005 and s 2m = (s^)^ ax = 0.010, according to the bounds of 
Eq. ([21]) ,p] In Fig. |2| we present the two branching ratios Brir — + cy) and Brir — > 3e) as 
functions of M 2 , for two fixed ratios M\jM 2 = 0.1 and 0.5, and for the CP phase £ = 0. 
We see from Fig. ^| that the LFV branching ratios in Model II are Br(r — > cy) ~ 10~ 8 and 
Br(r — > 3e) ~ 10 -9 . These values decrease relatively slowly when the parameters of the 
Dirac sector (a,b,c,d;l;) are moved away from the "optimal" values. This contrasts with 
Model I, where the maximal rates are reached only in a finely tuned region of parameter 
space. Our maximal values of Brir — > ej) agree with those of Ilakovac |3l|] - ours are by 
about factor three lower only because we took a different upper bound (s£ r ) 2 < 0.010 PT|). 

Fig. [3] shows the T-odd asymmetry At of Eq. (|20|) for the same choices of mass matrix 
parameters as in Fig. 0(b), but for £ = 7r/4 (solid line) and £ = 3tt/A (dashed line) (maximal 
CP violation). When all four heavy Majorana neutrinos are degenerate, there is no CP 
violation fl5| and At = 0. Also At = 0if£ = 0or7r/2. Notice that the maximal 5r's are 



reached for £ = 0, thus no CP violation; in the cases of Fig. | (£ = 7r/4, 3tt/A) the values of 
Br's are about one half of the corresponding ones in Fig. |2|(b) (£ = 0). The searches for 
maximal rates and for CP violation are in this sense complementary in their most optimistic 
cases. 

In the r — > \i processes, the maximal branching ratios in Model II are suppressed by an 
additional factor of (s^Lx/Oz )max = 0.002/0.005 w 0.4 [cf. Eq. (||)]. The approximate 
maximal branching ratios for the \x — > e LFV processes in Model II are obtained from the 
corresponding r — > e branching ratios by multiplying them with (s^ )max/( s L )Lx = 0-2 an d 
dividing by Br(r — > ^V^v T ) m 0.174. Thus, Br(fj, ej) ~ 10 -8 and Br(/i — > 3e) ~ 10 -9 , 



3 In contrast to Model I, we do not have the requirement of (s^ r ) 2 = (s 1 ^) 2 which followed there 
from the maximal {yt~ u T ) mixing condition (i = fi there), so that in Model II we can saturate each 
of the two upper bounds ( |2l| ) separately. 
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which is above the present experimental bounds 1.2 x 10 11 and 1.0 x 10 12 . Therefore, the 



maximizing conditions (|34[) -(35) cannot be met in this case. 



C. Expected numbers of events 



The explicit numbers of expected events in the considered processes depend on the way 
the r leptons are produced and on the luminosities involved. For example, the r pairs could 
be produced via e + e~ — > t + t~ close to the production threshold, or by sitting on a specific 
vector meson resonance V: e + e~ — > V — > t + t~. In this case, <r(e + e~ — > V) as a function of 
the CMS energy y/s can be approximated as a Breit-Wigner function 

,(.; e+e- _ V) = K^j-^-L—— , (36) 

where My and ry are the mass and the total decay width of the resonance. Constant K in 
Eq. (|36D can be fixed by invoking the narrow width approximation (nwa) formula 

a nwa (s; e + e- - V) = l2 ^^ V) 5{a - M v ) , (37) 

My 

where T ee (V) is the partial decay width for V — ► e + e~ . Integration of (|37j) over the variable s 
gives 12n 2 T ee (V)/M Vl fixing the constant K in ©: ^ = 37rr ee (\/)ry/M2 . The production 
cross section is maximal on the top of the resonance \fs = My : 

^e+e-^V^r+r-) 10 ** = a(e + e- ^ V^) max x w 12tt %^ ^^ t^ ■ ( 38 ) 

Multiplying this cross section by twice the branching ratio Br(r — > e^{eee)) we obtain 
the cross section for the process e + e~ — > V — > r + r _ — > e7(eee) + r. These branching 
ratios are ~ 10 _8 (10 -9 ) in Model II, as shown in Figs. [| For example, if the resonance 
is taken to be V = T(1S) (My = 9.46 GeV; T ee (V)/V v w r rT (V)/r y w 0.025), then 
cr(e + e~ — > V — >• r + r _ — > e7(eee) + r) would be about 2. (0.2) fb. For a luminosity of 
10 fb _1 /yr, this corresponds to 20 (2) events per year. Increased luminosities would give 
correspondingly larger numbers of events. 



V. CONCLUSIONS AND SUMMARY 



We investigated and compared heavy Majorana neutrino effects on lepton flavor violating 
(LFV) decay rates of tau leptons, in two popular models: (I) the interfamily seesaw-type 
model realized in the SM with right-handed neutrinos, and (II) the SM with left-handed 
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and right-handed neutral singlets. Further, we calculated a T-odd asymmetry At for r — > 
3/'. Model I is severely constrained in most of its parameter space by the actual eV-scale 
experimental upper bound on the light neutrino masses. It can give maximal LFV branching 
ratios Br(r — > /ry) ~ 10~ 9 and Brir — > 3fi) ~ 10~ 10 in a very restricted region of parameter 
space (where, incidentally, the CP- violating asymmetry At is zero), but otherwise it gives 
branching ratios many orders of magnitude smaller. On the other hand, in Model II the LFV 
branching ratios can be larger over a wide range of parameter values, Brir — > ej) ~ 10 -8 
and Brir — > 3e) ~ 10~ 9 , and At can reach values up to 5%. The results in Model II are 
insignificantly affected by the experimental bounds on the light neutrino masses. Model 
II can predict large enough LFV branching ratios (Br's) to be tested with near future 
experiments. The maximal LFV Br's in Model II are obtained when At = 0, and are reduced 
by about a factor of two when At ~5%. Model II thus indicates a certain complementarity 
in the searches of the LFV decay rates and of the associated T-violation. 
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Appendix A. APPROXIMATE MAXIMIZATION OF THE BRANCHING RATIOS 



The amplitude squared for the I — > I'-y LFV process is, according to Eq. 
mately proportional to 



), approxi- 



N L N L +N R 

|A| 2 oc|y>* J B M x/(0)+ £ B* lh B Vh xf{m 

t=l h=N L +l 



(A.l) 



where the first sum runs over the light, practically massless, intermediate neutrinos, and 
the second sum over the heavy neutrinos with masses (~ M\ ~ M 2 ). The function / 
depends of the mass of the exchanged neutrino; in the specific case, it is the loop function G 7 
appearing in Eq. (0). We now approximate the second sum by assuming that all the heavy 
neutrinos eigenmasses irth are the same: rrth = M. Then this, together with the unitarity of 
the matrix U (note: Bu = U* t ), implies 
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N L +N R 

\A\ 2 oc\ 

h=N L +l 



N L +N R 

U h2 U* hl \ 2 x\f(M 2 )-f(0)\^\ £ U h2 U* hl 

h=N L +l 



(A.2) 



Here we denoted the flavor index / of the heavier charged lepton as 2 and the index V of the 
lighter one as 1. 

The amplitude squared for the / — > 31' LFV process is, according to Eqs. an d 
(|T0|)-([l3|), more complicated. However, in the leading order in rriDrri^, and when there is no 
CP violation (when matrix U is real), it is straightforward to show that |v4| 2 is proportional 
to the same kind of combination ([A. 1|) . Thus, the proportionality ( |A.2| ) is approximately 
satisfied also for the I — > 3/' LFV process. 



The above proportionality can be supplemented by the Schwarz inequality 



a\ 2 <x\J2u h2 u* hl <J2\u h2 \ 2 xJ2\u hn \ 



(A.3) 



The equality is achieved only when there is a proportionality: 

_U h2 

T 

h=N L +l 



Urn 



u, 



111 



h=N L +2 



Urn 



h=N L +N R 



(A.4) 



Thus, the approximate maximal value of \A\ 2 , and thus of the LFV branching ratios, 
is achieved when the values of the heavy-to-light mixing parameters (s u L 2 ) 2 and (s^ 1 ) 2 are 
saturated according to the upper bounds fl2~l~|) and, at the same time, the mixing matrix U 
elements in the heavy-to-light sector satisfy the equalities (|A.4j ). 



In the seesaw Model I (with Nl = Nr = 2), the mixing matrix U elements in the 
heavy-to-light sector are Uh 2 = (m^m^h^ and Uhi = ( m M m DVi> where h' = h — 2. In 
our specific case of maximal mixing (a = c, b = —d, 8\ = 5 2 = tt/2) we have: U31 = a/ Mi, 
U 41 = -ib/M 2 , U 32 = -ia/Mx, U A2 = -b/M 2 ; the equality (gg ) is fulfilled; and (s^ 1 ) 2 = 



{a 2 /M 2 + b 2 /M 2 ). 



In Model II (with Nl = 2 and Nr = 4), with £ = 0, it can be shown, e.g. by using 
Mathematica, that the equality (|A.4|) is satisfied when ad = be. If in this case also the values 
of (s L 2 ) 2 = J2h \Uh2\ 2 and (s^ 1 ) 2 = J2h\Uhi\ 2 are saturated by the corresponding upper 
bounds of Eq. (|2l"D, then the approximate maximal branching LFV ratios are achieved. 
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FIG. 1. Maximal branching ratios for r — > ^7 (solid lines) and r — ► 3/x (dashed lines) as 
functions of M 2 in Model I, for a fixed ratio Mi/M 2 = 0.1 (a) and M 1 /M 2 = 0.5 (b). Mi and M 2 
are restricted to be above 100 GeV and below the perturbative unitarity bounds p^), indicated by 
the vertical line. The Dirac mass parameters are taken in the form (|3l])-(|32|) (and 5i = 5 2 = n/2) 
which give maximal branching ratios at any given M\ and M 2 . 
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FIG. 2. Maximal branching ratios for r — > e7 (solid line) and r — > 3e (dashed line) as functions 
of M 2 in Model II, for a fixed ratio M x /M 2 = 0.1 (a) and M x /M 2 = 0.5 (b) . Mi and M 2 are restricted 
to be above 100 GeV and below the perturbative unitarity bounds (p2|), indicated by the vertical 
line. The Dirac mass parameters a, b, c, d are taken in the form (p4|)-(|35|) which give approximately 
maximal branching ratios. Here we use the CP phase £ = 0. 
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FIG. 3. The T-asymmetry At in Model II, for the decay r — ► 3e as a function of M2, keeping 
Mi/M2 = 0.5 and adjusting the mass parameters in order to obtain maximal branching ratios [i.e. 
Fig. 2(b)]. The CP phase is taken to be £ = 7r/4 (solid line) and £ = 37r/4 (dashed line). Br(r — ► 3e) 
for £ = 7r/4 or 37r/4 is lower by about factor two in comparison to the case £ = of Fig. 2(b). 
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